Bose-Einstein Condensation in the presence of an artificial spin-orbit interaction 
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Bose-Einstein condensation in the presence of a synthetic spin-momentum interaction is consid- 
ered, focusing on the case where a Dirac or Rashba potential is generated via a tripod scheme. 
We found that the ground states can be either plane wave states or superpositions of them, each 
characterized by their unique density distributions. 
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The possibility of producing an artificial gauge field [l| 
opens up a new era in cold atom physics. The concept 
of gauge fields is ubiquitous in many branches of physics, 
from coupling of electromagnetic fields to charged parti- 
cles relevant to ordinary laboratory settings to fun- 
damental forces between elementary particles [3J. In ad- 
dition to the abelian gauge field already realized in [lj, 
there are numerous proposals to generate other artificial 
abelian and non-abelian gauge fields, with and without 
optical lattices, with the prediction of some exotic prop- 
erties [M3|. 

Two particular cases have captured much attention 
due to their relation to condensed matter physics: a 
Dirac-like term in the Hamiltonian and Rashba inter- 
action (3-[lo|. The former is relevant to systems such 



as graphene (15|, whereas the latter to surface states in 
topological insulators 16] and is also discussed frequently 
in the context of spintronics 



17 



and both of them in 
non-centrosymmetric superconductors f5~sj ] . 

While Dirac and Rashba interaction are often investi- 
gated in fermionic electronic systems mentioned above, 
in this paper we investigate the consequences of Bose- 
Einstein condensation (BEC) of Bosons in such artificial 
gauge fields. We investigate the possible ground states 
of this system, and show that one can have BEC into 
plane-wave states or states corresponding to a superpo- 
sition of plane waves, depending on the physical parame- 
ters. We further investigate the observable consequences, 
assuming that the gauge fields areproduced via the tri- 
pod scheme proposed in 0, [H, 0, HI- We shall see how 
the single particle terms due to the laser fields or ex- 
ternal potentials and the interparticle interactions deter- 
mine which ground state would be realized. Each ground 
state is associated with a characteristic density distribu- 
tion of the physical states involved in the tripod scheme, 
enabling one to distinguish these states experimentally. 

We mention here the closely related works of [rj [l3l[l4 |. 
[HI investigates BEC in the artificial gauge field in 
the experiment cited in [l4| considers also the 

Dirac/Rashba interaction as in the present work, but 
without reference to a particular production scheme and 
so did not discuss the potentials arising from the laser 
fields, nor the density distributions in the original atomic 
states. They also limit themselves to a less general form 



of the interparticle interaction included in this paper. [6j 
did not include the single-particle and interaction en- 
ergies to be discussed in detail below. They obtained 
a ground state qualitatively different from 13, and 
ours. 

For definitcness, we consider the tripod scheme pro- 
posed in Q. For convenience we shall first review it 
briefly. We consider a cloud which is either three- 
dimensional or confined by y-dependent optical poten- 
tials so that one only needs to consider motion along 
x and z. Three (almost) degenerate states |1 >, |2 >, 
1 3 > arc coupled to another state |0 > with light fields 
with strengths £2 1.2, 3. That is, we have the coupling term 
H a = -|0 > (Oi < 1| + fl 2 < 2| + 3 < 3|) + c.c. fl h2 
are chosen to be plane waves of equal strengths propa- 
gating in opposite directions, while involves a plane 
wave propagating perpendicular to the previous two with 
the same wavelength: we have Ox — -j^ls\x\9e~ tkaX , 

il 2 = -±ttsm9e +lk ° x , and fl a = flcos9e ik ° z , with Q and 
9 parameters specifying the coupling strengths. There 
are two linear combinations among the states |1 > — 13 > 
which are not affected by Hq. They are referred to as 
"dark states" which we shall choose to be 



|I>i > = 4= [I 1 > e+tk " X - |2 > e 



\D 2 > = -^cos<9 [|1 > e +lkox + \2> e - lkox ] t 
V2 



(1) 

■j ■//,-,:.;■ _[_ |2 \, ikox~\ c ikocos6 z 

sm9\3 > e -tto(i-c«9')* (2) 



Here we have generalized slightly [7j and include a yet 
undetermined parameter 9' in the z-dependent phase fac- 
tors in eq (J} and ([2]) . The value for 9' does not affect the 
actual physics, but we shall choose a special value for it 
for convenience below. The "bright" state \B >, orthog- 
onal to both 1-01,2 >, couples to |0 > with strength £2. 
We shall assume that the energy splitting of these two 
resulting states from |Z?i.2 > are large compared with all 
other energy scales considered below, hence |£>i,2 > are 
the only states physically relevant. 

Within the |-Di,2 > subspace, the single particle Hamil- 
tonian is given by 
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H ^i(^- 1 ) 2 + * + V < 3 > 

Here the matrix A has elements = i < D^VDv > 
(/i, v = 1 or 2) is the resulting (non-Abelian) gauge field, 
and $ M „ = ^ < Vigfl > • < B|V£>„ > 0. V in- 
cludes the terms that may arise if the internal energies 
or the potentials acting on the states |1 — 3 > are not 
identical. We shall provide more details below. A has 
components An — —koCosO'z, An — A21 — — k^cosOx, 
A22 = -k (cos6'-sm 2 6)z. With the choice cos6>' = s^£ ) 
we then have 



and J^- = 0. Obviously p y is zero and for simplicity we 
shall not write this component explicitly. The possible 
minima arc p = (0,±fco sl 2 ° ) = (0,±Pz), with energy 
E = Si (^), andp = (±k Q cos6,0) = (±p°,0), with 

energy E = — ■^cos 2 9. Hence if cos9 < \/2 — 1, the 
minima are at p — ±p®z, whereas if cos# > \f2 — 1, 
the minima are at p = ±p%x. If cos# = y/2 — 1, we 
have the very special case that all momenta given by 
{p 2 x +P1) 1 / 2 — k Q (y/2 — 1) are degenerate. We shall not 
deal with this very special circumstance in the present 
paper. 

If we have condensation in one of these four minima, 
the corresponding wavefunctions in |-Di,2 > space are 



-knCosOa,, 



— Sin Wer, 
2 



(4) 



where <J x ,y,z are Pauli (pseudospin) matrices acting 
within the \D\ >,|i?2 > space. With this A, we see 
that we have an anisotropic Dirac-like term in the Hamil- 
tonian ^ (cos6(r x p x + sm 2 8 a z p z ^j , with the special case 

cosd = \/2 — 1 = where this term becomes isotropic 
in the x-z plane 0| ■ We shall however not assume such a 
special value for 9 below. 

In the above vector potential A, x and z result 
from the x and z dependence of VL\ t 2 and re- 
spectively. One can thus produce a Rashba-like term 

^ (cosda x Pz — ^^^Y^VzPxj if one replaces x by z in Q,^ 
and z by —x in SI3, that is, changing the directions of the 
lasers or relabeling the coordinates. Alternatively, one 
can also use a different choice for \Di^ > corresponding 
to a "spin-rotation" Q. For defmiteness, we shall con- 
tinue to deal with the Dirac-like term in our Hamiltonian. 

$ can be easily found to be 
^sin 2 ^ ( 1+c 2 os2e + ■ The first term is a 

momentum independent scalar which we shall drop. 
In the presence of a potential (or internal ener- 
gies) Vi,2,3 on the states |1 — 3 >, V is given by, 
apart from a constant which we shall again drop, 



V = 



{Vi - V 2 )cos9a x + {V! +V 2 - 2V 3 ) s -^a 



and V act like Zeeman fields in the space \D1_2 >■ 

We shall first consider BEC in the presence of the 
gauge potential A, but ignoring V, and interparticle 
interactions for the moment. These would be included 
later. The Hamiltonian is then simply the kinetic energy 

is sim- 



term K = 7r—(p 

2m V* 



A)2=-L,?-lAf 
ply given by the constant k$ ^cos 2 # ^ 



2m 



which we 



shall drop below. K can easily be diagonalized giving the 



energies E = §^ 



± 



nl/2 



PI 



BEC 



should occur in the state of the lowest energy, hence we 



make take the negative sign and find the p where M^- = 



* 4 (r) 




1 
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3_ ikocosOx 

V2 \ i 



V2 



4_ — ikocosOx 



1 



(5) 
(6) 
(7) 
(8) 



where $1-4 are complex numbers. We shall call them 
states 1-4 (not be be confused with those which enter the 
tripod scheme). We note here that even though these 
are plane wave states, they carry no current since we 
have demanded d ® E = 0. Alternatively, we should note 

that the velocity operator v is given by hence v = 

. One can easily 



~(p-A) = £ + cos6a x x 

m^r ' m m \ x 



<T Z Z 



verify that the expectation values of the velocity is zero 
in each of the plane wave states ([5][5]) above. 

Since in each case we still have two degenerate min- 
ima, we must consider other terms in the Hamiltonian 
to determine which BEC would occur. We first consider 
the other single particle terms in the Hamiltonian, that is 
V+$ = H/j = —h-a, assuming for the moment that they 
are dominant over the interparticle interactions. Due to 
the above discussion, we shall consider the case where H/j 
has only a XlZ terms, With finite h, the relevant branch of 
the spectrum is 



EL 
2m 



E=r 



cos0p x 



sin 



21 V2 



(9) 

We now need to find the value of p where E is a min- 
imum. In general this is very complicated for general 
h. In the limit of first order in h, one can easily ver- 
ified that the energies of the above four states become 
E = ± h z and E = ~^cos 2 6 ± h x . Hence 



if cosO < y/2 — 1 and h z < (>)0, then state 1 (2) should 
be realized, whereas if cos(9 > y2 — 1 and h x < (>)0, 
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states 3 (4) should be realized. The above statements as- 
sumed that h is sufficiently small so that the ordering of 
energies are not changed except the lifting of degenera- 
cies. Also, strictly speaking, one also needs to consider 
the modification of the wavefunctions in eq ([S][H]) due to 
h, as both the momentum and the "spin" wavefunction 
at which the minimum occurs are modified. We however 
would not give these rather lengthy formulas here. (See 
however near the end of this paper below). 

Let us at this point give some physical properties re- 
lated to these states, ignoring for simplicity the modifica- 
tions of wavefunctions just mentioned. If cos# < \/2 — 1 
and h z < 0, then (0 is the wavefunction in \Di 2 > space. 
Using ([1} and ©, one easily finds the wavefunctions in 
the original |1 — 3 > basis: 

^i,2(f) = ^cos6e ±lkox e tkosin2ez (10) 
v2 

^ 3 (f) = -$ lS in(9e- lfeoCOs2ez (11) 
and so the corresponding particle densities are 



cos 2 



|V 3 (r)| 2 = |$i| 2 sin 2 



2 



(12) 
(13) 



Hence IV'iWI 2 = I ^(r)] 2 (again ignoring corrections due 
to h, a statement which we shall not repeat). 

If cos# < \J~2 — 1 and h z > 0, the system condenses 
into state 2 (eq J6j) , we have instead 

2 ±ikox 



^ 2 (r) = ±^=e="'" 



(14) 

= (15) 
and with the corresponding particle densities 



l^i, 2 (r)| 5 



1$, 



\M^)\ 2 = o 



(16) 
(17) 



which is very different from the state 1. 

For cos0 > \/2 — 1 and h x < we have condensation 
into state 3, where 



</, 12 (f) = -^[= F l + cos^e ± *°( 1±cos ^ a: e ifeoM ^ z (18) 
v2 



The densities are 

Wl,2W = l $ 3| 



2 _ I* ^^TcosO) 2 



2 sin 2 6» 



\Mr)\ 2 = I* 



(20) 
(21) 



For cos# > \/2 — 1 and h x > we have condensation 
into state 4. We have 

fa(f) = _^ s i n 6l e -ifeDcosfe e -ifeo(l-^) 2 ( 23 ) 



with densities 



|^i,2(r)| 2 = |* 4 | S 
l^ 3 (r)| 2 - |$ 4 | : 



(1 ± cos60 2 



, sin 



(24) 
(25) 



which are those of state 3 with |1 > and |2 > inter- 
changed. In passing, we remark here also that, as indi- 
cated by the wavefunctions ^1,2,3 (?) above, the atomic 
states |1 >, 1 2 >, |3 > are each associated with their 
characteristic wavevector, which can also be measured 
by time-of-flight experiments as in [l|. 

The above is for when the single particle terms dom- 
inate and determine the ground state. Generally, if the 
energy differences between the above plane wave states 
are sufficiently small, we can have condensation into su- 
perposition of plane wave states. The energetics is anal- 
ogous to that of superfluid mixtures [l9} as we shall see. 
To be specific, let us consider the case where we have 
cos# < y/2— 1, with energies of states 1 and 2 sufficiently 
lower than 3 and 4 so that the later two need not be con- 
sidered. We are left with considering condensation into 
states described by 1 and 2 in eq (p} and (J2j) . (We are 
now considering normalized single particle states, that 
is, without the $1,2 coefficient). Consider the scattering 
between two particles, each one in the state 1. Since the 
total momentum is 2p , (here po = (0, ko^j^-)), the only 
relevant term in the low energy sector must correspond to 
that that the outgoing particles are also in state 1. Simi- 
larly, if there is one particle from each state 1 and 2, then 
the total momentum is zero and the relevant low energy 
process also correspond to outgoing particles with one 
each in state 1 and 2. Let dp .i and d-p 0:2 be the opera- 
tors corresponding to the (normalized) plane wave states 
in states 1 and 2. The general form of the Hamiltonian 
in the low energy subspace is 55n^ ,i d ]? 0i Wpo,i rf Po,i + 

9l2<l i ; ol d}_ i!a2 d^ .2dp ^ + |522rfLp- 0i2 d -po,2 rf -Po,2 rf -po,2, 

where c/ij are coefficients. Other terms in the Hamilto- 
nian are off-shell and can be dropped when we eventu- 
ally take the mean-field approximation below, gtj can 
be evaluated once the interaction Hamiltonian in terms 
of the states 1 1 — 3 > in the tripod scheme is known. 
There are two many possibilities so we do not provide 
the formulas here. We simply remark here that if the 
interaction among the particles |1 — 3 > are all identical, 
then gii = g 12 = g 22 

Within the mean-field approximation, we replace the 
operators ^±^,,1(2) by the amplitudes $1,2, obtaining the 
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interaction energy E mt = ±<7 n |$i| 4 + gi 2 |$i| 2 |$ 2 | 2 + 
5ff22| ( i ) 2| 4 - Also, in the presence of the bias field Hh, 
we have energy — h z (\$i\ 2 — l^l 2 ) to linear order 
in h corresponding to the discussion below eq ©. Note 
that there is no such term such as — h x (<&* 1 § 2 + $2^1) 
There cannot be a term in the form of —h x (dg o x d-p~ .2 + 
dt ,<!-»*„ 1 ) due to momentum conservation. 

Including also a chemical potential fj, gives us then the 
total energy E tot = -(//- h z )\$i | 2 - (fi+ h z )\<& 2 \ 2 + E int 
as in the case of a mixture of two species 1 and 2, with 
effective chemical potential for component 1 (2) being 
Mi, 2 = /iT h z - Note that therefore there are no terms 
which would depend on the phase difference between 
<&i 2 ■ We shall have another point of view of this be- 
low. The phase diagram of a two component mixture 



is well-known [19(, and we shall not repeat those results 
here. The state corresponding to a "mixture" between 1 
and 2 has a wavefunction in the \D± t 2 > space as a the 
superposition of {2} and ([6]), that is, 

and likewise for their projection into the states 1 1 — 3 > . 
We can understand easily why the energy is independent 
of the relative phase between "Ji^- Changing the relative 
phase between these two complex numbers can simply be 
reabsorbed by a shift of the origin for z. Below, we shall 
for simplicity assume that both $12 are real. 
The densities for this "mixture" of 1 and 2 are 



l^i,2(r)| 2 

1^(01 



1 r 2 

2 [ 



|$i| 2 + I $ 2 1 2 ± 2$i$ 2 cos6»cos(£:osin 2 6»,2)] 



2 = ISilW 



(27) 
(28) 



The coefficients of |$i,2| 2 are the same as eq ([T2"j) - ([TT]) . An interference term oc $i$2 arises from the wavevector 



difference fcosin 92 between the $1^2 components in eq (|26p . Similar discussion holds for a 
4. The densities are 



mixture" of states 3 and 



liM*Dl a 
l^)| 2 



- cos6») 2 |$ 3 | 2 + (1 + cos6i) 2 |$4| 2 - 2<J>3$ 4 sin 2 l9cos(2fc cos6>20] 
[|$ 3 | 2 + |$ 4 |2 + 2$3$ 4 cos(2fc cos(9x)] 

I 



(29) 
(30) 



The densities thus in general oscillate in space, as has 
been discussed in 13l 14|. 

Next we would like to consider in more detail the case 
of H/j = —h z a z , without restricting to small h z . Recall 
that 3? contains such a component even in the special case 
that V = 0, and also that h z < in that case. Due to the 
lack of space we would mainly just state the results. First 
we ignore interparticle interactions. If cos9 < s/2 — 1, the 

) is the energy minima 

For 



state 1 (2) with p = (0, ±fco 
if h z < (>)0 with energy E 

COS9 > y/2 

with D{9) e 
range for 9) 



sin 2 9-\ 



COS 



< h c where h c = %D(9) 
(note D > in the above 



1, and if \h z \ 
, then the energy minima occurs at (±p~.,p~) 



where = fc cos6»[I- (%-) 2 ] 1/2 andp° = -h z ?fJL. Here 
h z = hz/ik^/m) is a dimensionless measure of the energy 
h z . Note there are two degenerate minima of opposite 
p x , and sgn(p° z ) — —sgn(h z ). We shall call these states 

k 2 

5 and 6. The energy is given by E — — ^ 



cos V + 



At Th c , p becomes (0, ±k ^j^) and thus merge with the 
other minima stated before. For \h z \ > h r , the minima 



is at p = (Ojifco 5 ^) accor ding to h z < (>)0. 

It is convenient to define the quantity /3 = 



tan 



-h z /D 



l-(h z /D) 2 



13 = if h z = and p = -sgn(h z )Z 

at \h z \ = h c . The wavefunction in \D\^ > space for state 
5, with wavevector (p x ,p° z ), is 



-(cos§ — sin |) 



cosf 



sinf 



(31) 



whereas that for state 6 with wavevector (— p%,p° z ) is 



v2 



(32) 



These two states are orthogonal due to their different 
wavevectors (±p x ,p®), but their "spinor" part has finite 
(sin/3) overlap. At h z = 0, = 0, p x = k Q cos9, p° z = 0, 
states 5 and 6 are identical respectively with state 3 and 
4 . At h z = —h c , P = § j and as stated, p x = 0, p\ = 
k ^—, they become identical with each other and with 
state 1. 

The densities corresponding to state 5 are 
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- [l + cos 2 # - sin/3sin 2 6> =F 2cos#cos/3] 

; 

-sin 2 6»(l + sin/3) 



(33) 
(34) 



The densities of state 6 are those of 5 with states |1 > 
and 1 2 > interchanged. 

For a superposition between states 5 and 6, the densi- 
ties have three contributions, the terms proportional to 
l^s^l 2 are the same as those in eq (|3"3"j) and (jM]) . In 
addition, there are the interference terms 



the Aspen Center for Physics where this study was mo- 
tivated. 



[sin 2 (9 - sin/3(l + cos 2 (9)] cos(2p°x) 



for 1-01,2 (^0 1 2 , and 

$ 5 $ 6 sin 2 6>(l + sin/3)cos(2p°a;) 

for |?/>3(r)| 2 . Note that they reduce to the appropriate 
limits eq (|29]) flU if /3 = 0. At /3 = §, we get eq (JT2J) 
and ([13]) with $i -> <5> 5 + $ 6 . 



If V = 0, we have h z = 



Hence \h z \ = h r 



at cos6< = m 0.518, corresponding to w 0.326?r. 

Hence, due to this finite h z , we have states 5, 6 or their 
mixtures if 9 < 0.3267T, and state 1 if otherwise (assum- 
ing that interaction does not change the ordering of the 
energies). 

In [ij], standing waves for f2i :2 are considered instead 
of plane waves. This situation can be treated in a similar 
manner as in this paper. The corresponding results can 
be obtained by a simple unitary transformation among 
the atomic states |1 > and |2 >. 

In conclusion, we have considered Bose-Einstein con- 
densation in an artificial spin-orbit field. In general BEC 
occurs in a state with finite wavevector, or their super- 
positions. Due to the momentum-pseudospin coupling, 
each state is characterized by their unique density distri- 
butions. 
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